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Abstract 

We investigate a family of lattice models with manifest M=2 supersymmetry. The models 
describe fermions on a ID lattice, subject to the constraint that no more than k consecutive 
lattice sites may be occupied. We discuss the special properties arising from the supersym- 
metry, and present Bethe ansatz solutions of the simplest models. We display the connections 
of the k = 1 model with the spin-i antiferromagnetic XXZ chain at A = —1/2, and the 
k = 2 model with both the su(2|l)-symmetric tJ model in the ferromagnetic regime and the 
integrable spin-1 XXZ chain at A = —1/^/2. We argue that these models include critical 
points described by the superconformal minimal models. 



1 Introduction: supersymmetric lattice models 

In the analysis of quantum mechanics and quantum field theory, one is often guided by the 
presence of symmetries. Among these symmetries, supersymmetry stands out as particularly 
powerful: supersymmetric theories enjoy special 'magical' properties and in many cases super- 
symmetry leads to considerable computational simplification. 

Supersymmetric field theories in two (2-1-0 or H-l) dimensions are particularly interesting. 
They have manifold applications, ranging from their use as world sheet theories in (super-)string 
theory to their role as continuum limits of critical models of 2D statistical mechanics or ID 
quantum chains. For example, the field theory describing the (experimentally realizable PP) 
antiferromagnetic XXZ spin chain at anisotropy A = —1/2 has M={2,2) supersymmetry. 

Other lattice models, such as the two-dimensional tricritical Ising model along the first-order 
transition line, are described by supersymmetric field theories in the continuum limit as well. 
In none of these cases is it obvious how the supersymmetry manifests itself in the underlying 
lattice model. To understand this issue, and to show how supersymmetry improves our under- 
standing of lattice models directly, we introduced lattice models with explicit supersymmetry ^ . 
These models describe hard-core fermions on a general lattice or graph, and the supersymmetry 
generators are related to operators that create or annihilate these fermions. We showed that the 
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simplest of these models defined on a ID chain is critical, and is closely related to the XXZ chain 
at A = -1/2. 

In the present work, we greatly expand upon the work of [3j|. We introduce a series of 
supersymmetric lattice models, the simplest one that of Here we we limit our study to one- 
dimensional (open or closed) chains, but the models (if not the results) can be easily generalized 
to higher dimensions. Analyzing the models in detail, we pay particular attention to the following 
aspects: (i) special features, such as zero-energy ground states, implied by the super symmetry, 
(ii) integrability by Bethe ansatz, (iii) mappings to other models, and (iv) continuum limit and 
relation with models of AA=(2,2) superconformal field theory. In section 2, we define a series of 
models (denoted M^), and compute their Witten index. In section 3, we discuss the simplest 
model Ml in detail, extending the results of [Hj- In section 4, we show that the model M2 at 
a particular coupling is closely related to the stt(2|l)-supersymmetric tJ model, thus providing 
an intriguing relation between spacetime supersymmetry and global symmetries with fermionic 
generators (which are often called supersymmetries in the condensed- matter literature). We also 
show that at a different coupling, M2 is related to the spin-1 XXZ model with couplings tuned to 
be integrable. The latter model and the model Mi are described in the continuum limit by the 
first two minimal models of superconformal field theory. In section 5, we argue that a particular 
coupling of the model is described by the kth superconformal minimal model. 



2 A family of J\f=2 supersymmetric lattice models 

2.1 The models M^.: general structure 

To define the models M^,, we consider fermions on a ID chain, subject to the condition that 
at the most k consecutive lattice sites may be occupied. This defines a Hilbert space Tik, with 
subspaces defined by the condition that the eigenvalue of the fermion-number operator F 
be equal to /. The supersymmetry generators Q = Q~ and = act on the constrained 
Hilbert spaces TCk as 

Q-:ni^^^nr\ Q-:n^P^nr'' (i) 

so that 

[F,g±]=±Q±. (2) 

They are nilpotent: 

(Q+)2 = (Q-)2 = . (3) 

One of the most important characteristics of models with spacetime supersymmetry is that the 
algebra of the supersymmetry generators involves the Hamiltonian H. For J\f=2 supersymmetry, 
the relation is 

H = {Q+,Q-]. (4) 

This definition, together with the nilpotency of the supersymmetry generators, implies that H 
commutes with both the fermion number and the supersymmetry generators 

[F,F]=0, [Q+,H] = [Q-,H]=Q . (5) 

The algebraic structure summarized in the equations ((21), (jSl) ® and © is characteristic of 
supersymmetric quantum mechanics or J\f=2 supersymmetric quantum field theory. 
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To further specify the models Mk, we start by introducing creation and annihilation operators 
c|, Ci for unrestricted fermions. They obey the usual anticommutation relations 

{q,c]} = 6ij {ci,Cj} = {c|,ct} = 0, 

so there is at most one fermion per site. The supersymmetry generators are particular combina- 
tions of the fermion annihilation and creation operators. There is some freedom in defining the 
on Tij. while maintaining the above algebraic structure. This freedom will lead to adjustable 
parameters xi, . . . , Xk-i in the definitions of the models M^. 

For k = 1, adjacent fermions are forbidden; i.e. the Hilbert space is constrained so that 
nearest-neighbor sites on the lattice cannot be occupied. With this, the fermions can be thought 
of as hard-core dimers centered on site i, with the dimer covering the sites i — ^ and i + ^ on 
the dual lattice. Introducing the projection operator Vi = 1 — c\ci, we can write di = 'Pi-ic\Vi+i 
and define the supersymmetry generators by 

N N 
1=1 i=l 

The resulting Hamiltonian will be displayed, and its eigenstates found, in section |21 

Putting A; = 2, we can define d\ = Vi-ic\Vi+i, e\ = Vi-2cl_iCi-icl'Pi+i + 'Pi-ic\c\^^Ci+i'Pi+2 
and define supersymmetry as 

Af 
i=l 

The condition (Q^)^ = is satisfied for all yi, y2, leaving the ratio xi = yi/y2 as a non-trivial 
adjustable parameter. 

For general fc, one may extend these results to construct a family of nilpotent supersymmetry 
operators with adjustable parameters xi,...,Xfc_i. This construction goes as follows. Define 
operators d[a,6],j) d] w ^ with o, 6 = 1, . . . , /c, 6 < a, and with d[a,fe],i taking out a fermion at a site 
i, which is at the 6-th position of a length-a string of consecutive occupied sites. The definition 
of c?[a,6],i includes the appropriate minus sign for bringing the annihilation operator to site i. 
Weighing the contribution of d[a,fe],i to the supersymmetry generator with ya^b, one arrives at 

= ^^ya,bd[a,bli ■ (6) 

i a,b 

Imposing then that (Q")^ = leads to recursion relations on the coefficients ya,b, which can be 
solved in a closed form that is reminiscent of Newton's binomium. Writing ya = ya,i, the solution 
reads 

yaya-i ■ ■ ■ ya-b+l ,r,^ 

ya,b = • (7) 

yb-.iyb-2 • • • 2/1 

This leaves yi,y2, ■ ■ ■ ,yk as free, adjustable parameters in the definition of a nilpotent supersym- 
metry operator on the space TCk- The resulting Hamiltonian depends non-trivially on all /c — 1 
ratios Xi = yi/yk, i = 1,2, . . . ,k - 1. 
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2.2 Witten index and ground-state degeneracy 

The algebraic structure @ and © imphes the fohowing simple properties of the spectrum 
ofH HI: 

• all energies satisfy E > 0, 

• all states with E > can be organized into doublets of the supersymmetry algebra, 

• states are annihilated by both and if and only if they are zero-energy eigenstates 
of H. 

The Witten index of a supersymmetric theory is 

Wk = Ttn,{-lfeM-PH). (8) 

Because all states with E > appear in doublets with the same energy and fermion number 
differing by 1, their contribution to Wk cancels. Thus Wk counts the number of ground states 
of the theory, weighted by (—1)^. A non-zero Witten index requires the existence of at least 
\Wk\ ground states with E = 0. Note that Wk is independent of /3, and in fact H itself, as long 
as supersymmetry is preserved. The number of ground states in the theory can change as the 
parameters change, but Wk cannot. This means ground states can only appear or disappear in 
pairs with opposite values of (—1)^. 

Since Wk is independent of /?, we can evaluate it at /? ^ 0, so that all states contribute with 
a factor (—1)^. We compute this by employing a recursion matrix to pass from to iV-|- 1 sites, 
while keeping track of the configuration near one open end. The characteristic equation for the 
matrix pertaining to the space TCk is 

A^+i_A'= + ...-(-l)^ = , 

which implies that A'^"^'^ = (— 1)'^. This means that the Witten index Wk{N) for an A^-site chain 
obeys ^^^(A^ + A: + 2) = {-l)''Wk{N). We obtain (for N = l,...k + 2) 

• 7ik on an A^-site open chain, 

, , f forl<N<k 

Wk(N) = { ^ - - (9) 

^ ^ ioT N = k + l,k + 2 ^ ^ 

• for TCk on an A^-site closed chain with periodic boundary conditions. 



(-1)^+^ forl<A^<A;-M 
(-l)*=(A:-h 1) foTN = k + 2 



Wk{N) = {) ' ^ ,r , (10) 



For the case with periodic boundary conditions, one can refine this result by restricting the 
index to sectors with a specific eigenvalue t of the shift operator T. As = 1, these eigenvalues 
are of the form t = exp(27ri//Af), / = 1, 2, . . . , A^. For the case N = n{k + 2), where \W\ = k + l 
we find 

W = for t = (-1)^+1 expr^j") with j = 1, . . . , A. + 1 (11) 
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and W = otherwise. For all other A'^ 

W = (-1)^+1 at t = (-1)^+1 , (12) 

and W = otherwise. 

These values of the index are independent of the Xi: as noted above, Witten's index is inde- 
pendent of all supersymmetry-preserving deformations of the theory. We will show in subsequent 
sections that in two simple cases, Mi and M2[x = 0], the eigenstates of H can be found using the 
Bethe ansatz. In these two cases, and for M2[x = V2] we also have non-trivial mappings relating 
the supersymmetric chain to other integrable models. The models Mk[xi . . .] do not appear to 
be integrable in general, however. 



3 Mi: a supersymmetric model of hard-core fermions 

3.1 Hamiltonian 

The Hamiltonian for the k = 1 model is found by working out the general expression eq. @ for 
the supersymmetry generators defined in eq. (|2.1j) . Denoting di = Vi-iCiPi+i to be the hard-core 
fermion, we have 

TV 

Hi = Y, + dl+A + Ad^+i. (13) 

2 = 1 

Henceforth we take periodic boundary conditions, so the index i is defined mod N . The potential 
counts the number of "empty pairs", vacancies two sites apart. Adding a fermion more than two 
sites from any other fermions decreases the potential by 2, since there are now two less empty 
pairs. It is convenient to rewrite the Hamiltonian as 

N 
i=l 

In this formulation, the Hamiltonian has a constant term, a chemical potential of 2 per fermion, 
a hopping term, and a repulsive potential for fermions two sites apart. This looks like a lattice 
version of the Thirring model with a hard-core repulsion, and a next-nearest-neighbor repulsive 
potential. 

We remark that the Hamiltonian H14|) has been introduced and analyzed in the literature; it 
is a special case of the so-called constrained XXZ model discussed in 0. These papers do not 
mention the supersymmetry properties which form the focus of our analysis here. 

3.2 The correspondence with the XXZ model 

The Hamiltonian p4jl bears a simple relation with the XXZ Hamiltonian with A = —1/2 



Here the a are Pauli spin matrices. The transformation can be effected on a periodic chain as 
follows. Each edge between empty sites is replaced by an up-spin, and each occupied site with its 
adjacent edges is replaced by a down-spin. The length L of the XXZ chain is therefore related to 
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the length of the fermion lattice as L = N — f. To transform the Hamiltonian, observe that the 
hopping term simply allows a fermion to move to a neighboring site, provided its other neighbor is 
empty. The effect is the same as that of the terms l/2(crj cr^_(_ +crjc7j^]^) in the XXZ Hamiltonian. 
The diagonal term in H14() counts the number of empty second-neighbor pairs and translates into 
the operator ~ "^P/^ + (1 + + ^j+i)/'^' which counts the number of down-spins and 

the number of nearest neighbor up-spin pairs. On a periodic chain this expression is identical to 
(3L -|- <7|o"|_,_]^)/4. This gives, up to an additive constant, a map between the Hamiltonians, 
Hi ^ Hxxz-i/2 + 3L/4. Note, however, that the constant 3-L/4 in the spin chain would not be 
constant in the fermion model, as L = N — f . 

In spite of this reasoning these Hamiltonians do not have entirely the same spectrum. One 
reason for this is that the size of the Hilbert space is not the same between the two models, even 
in the corresponding sector. The size of the sector (A^, /) in the fermion chain is (for < / < N/2 
and N/f not an integer) 

' N-f-1 
/-I 



N 

while the size of the corresponding sector in the XXZ chain is 



N-f 
f 

The second difference is the fact that if a fermion on position N hops to the right to position 1, 
the weight of the resulting state will have an additional minus sign if the total number of fermions 
is even. 

To relate the spectra of the two models we first identify the missing states of the fermion chain. 
The rules given above give a one-to-one correspondence with states of the XXZ chain, for all states 
in the fermion chain in which the leftmost position is empty. The action of the Hamiltonian ()14() 
within this subspace is well represented by the expression (|15|) . but of course the Hamiltonian 
will also hop particles onto the leftmost site. In order to adjust the Hamiltonian (|15|) to account 
for jumps outside the subspace, we have to restrict it to sectors where the translation operator 
T has eigenvalue t. Within this sector each state can be identified with t {t~^) times its left 
(right) translated version. When in the fermion chain the leftmost position is occupied by a 
left hop, we identify this state with t times its left translate. With this identification, the hop 
corresponds precisely with a spin exchange between the leftmost and the rightmost position, 
leaving the rest of the configuration unaffected. If the number of fermions is even an additional 
minus sign should be taken into account for the cyclic permutation of the fermions in this move. 
The result can be emulated in the in the XXZ Hamiltonian by a twist in the periodic boundary 
conditions: cr^^-^^ = {—l)-l'~^^taf and cr'[_^^ = {—l)-l'~^^t'^^a^ . Notice that with this boundary 
condition the eigenvalues of T in the spin model are not the usual roots of unity = 1, but 
those of the corresponding fermion chain = 1. With this modification the spectrum of the XXZ 
Hamiltonian in the appropriate momentum sector corresponds precisely to that of the fermionic 
chain. We will confirm this result below by showing that the Bethe-ansatz equations of the 
fermionic chain of sites are the same as the XXZ chain of length L = N — f \i one includes a 
twist factor t{—l)f~^^ in the latter. 
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3.3 The Bethe ansatz 
3.3.1 The Bethe equations 

The Bethe ansatz computation follows closely the original computation of Bethe for the Heisen- 
berg spin-chain. For a review, see e.g. [SI EI. An eigenstate with / hard-core fermions is of the 
form 

= ^ . . . v)« . . . 4^10) (16) 

where we order 1 < ii < i2 — 1 < is — 2 . . . . The state |0) is the vacuum state with no fermions, 
not the ground state. First we construct eigenstates of the translation operator T, which sends 
Tdt (it • • . = (it (it If there is no fermion at site A^, then this can be done simply by 

taking 

ip{ii,i2,...if) ~ ^1^2' •••/^/ • 

but this prescription does not work if = N. The way to fix this is to sum over permutations 
of the Hi. Denoting P to be the permutation (PI, P2, . . . , Pf) of (1, 2, . . . /), we define 

(p{ii,i2,...if) = ^App}^^ii%...ii^j. (17) 
p 

The periodic boundary conditions give a constraint on the amplitudes Ap. For to be an 
eigenstate of T, the amplitudes must be cyclically related as 

^P/,P1,P2,...P(/-1) = ^p/(-l)'^~"^^Pi,P2,...p/- (18) 
The eigenvalue t of T is then given by 

/ 

t = J{{^^.r' (19) 

i=l 

For periodic boundary conditions, untwisted in the fermionic basis, we set = 1. Defining the 
bare momentum pi as = exp(ipj), the relation (|18|) corresponds to quantizing the momentum 
of a particle in a box. 

Equation (|17)) is Bethe's ansatz for the eigenstates. There is no a priori guarantee that it will 
work, but in this case it does. The /ij and the amplitudes Ap are found by demanding that (j)^^^ 
be an eigenstate. Let us first study the case with two fermions. We have 

(^(2)= ^ V'(^l,^2)« 

«1<«2 — 1 

where the sum on ii and i2 is restricted to be 1 < ii < ^2 — 1 ^ A — 1. The Bethe ansatz is that 

92(n,«2) = ^vlpMpi^P2 
P 

Requiring (p^"^^ be an eigenstate of T means that 

A12 = -^21/if , 
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with t = (/Ui/i2) ^- Operating with the Hamiltonian 1)141) on one finds 

N 
i=l 

^ v'Cn, .2) + dl-A, + + • (20) 



«l<i2 — 1 



Shifting the latter sums by ii — > ii — 1, ii ^ ii + 1, 12 ^ ^2 — 1 and ^2 — > ^2 + l, respectively, 
yields 



N 
i=l 

where 

E = N-2f + ni+iJ2 + (/ui)"^ + (/i2)"^ 

The other terms must vanish for (p^^^ to be an eigenstate of H. They involve fermions two sites 
apart, and are 

Xi = ip{i, i + 2)- ip{i, i + l)-ip{i + l,i + 2). 

The first contribution to Xi comes from the potential energy. The other two are a little bit more 
subtle. Because of the hard-core repulsion, there are no terms in (|2fl|) of the form (p{i,i + l)dldl_^^. 
After the shifts in ii and 12, such terms would result in necessary contributions to Ecji^^'^ . However, 
since they are lacking, we need to add them into Ecf)^^^ and therefore subtract them off from the 
X^. 

The terms Xi must all vanish if (j)'^^^ is to be an eigenstate: for all i 

Xj = ^ tIpiHp2 (/ip2 ~ ^P2 ~ MP1^P2) = 

p 

Each of these vanishes if 

A2I _ fJ-2 fJ-2fJ'l + 1 - /i2 

^12 fJ-i /^2/^i + 1 - ^1 ' 
Combining this with the earlier condition yields 

Av2 _ _ N _ /^l(^l/"2 + 1-^1) 
A2I /U2(Ail/U2 + 1 - ^2) ' 

Using the fact that t^^ = 111112 gives the Bethe equation 

which holds for {k,j) = (1, 2) and (2, 1). One solves these two equations for fii and ^2 subject to 
the constraint t"^ = fiifj,2- Then the corresponding eigenstate is found (up to an overall constant) 
by substituting the values of fii into the equation for A2\l A\2. In general, one must also restrict 
A^i / /^2, or else 4>^'^^ vanishes. There is one exception to this: if //i = /U2 = exp(ibi7r/3), then the 
Xi vanish without imposing a condition on the Ap. These sort of states are discussed in detail 
in where it is explained how the Bethe ansatz gives a complete set of states. For example, 
one can easily check that ix\ = [12 = exp(ibi7r/3) yields eigenstates of energy E = N — 2 when 
is an odd multiple of 3. 
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This computation can be generalized to all /. The constraint that (/'^■^^ be an eigenstate is 
basically the same as the two-fermion case: one requires that 

ip{ii,i2, ...i,i + 2,... ,if)-ip{ii,i2, . . . + . . . ,if) - ip{ii,i2, . . . ,i + l,i + 2, . . . ,if) = (21) 

for any choice of ii,i2, ■ ■ ■ if and i. The trick to make this vanish is to consider the permutation 
P', which differs from P only in that Pk and P(k + 1) are reversed, i.e. P' = P1P2 . . . P{k — 
l)P{k + l)PkP{k + 2)...Pf. Then ^ can be written as 



— i^pii^P2 ■ ■ ■ i^pki^p(k+i) ■ ■ ■ ^^pf (/^p(fc+i) ^p(fc+i) f^pkf^p{k+i) 



This vanishes if for all P and k 



with 



At 



Ap 

g{a,b) = 



gifJ-Pik+i)! fJ-Pk) (22) 
a(ab + 1 — a) 



b{ab + 1-6) 

One can think of g as the bare S-matrix describing the phase shift when two fermions are inter- 
changed. The fact that g{a,a) = — 1 means that the wavefunction vanishes if two of the fij are 
identical (except for fij = exp(ibi7r/3), as discussed above). 

To find the /Uj, we impose the boundary condition H18|) . Note that 



g{fJ-Pf,HPi] 



>^P/,P1,P2,...P(/-1) ^ ^Pl,P/,P2,..P(/-l) 
Api^P2,...Pf Api^p2...Pf 

f 

= Y{9{l^Pf,f^Pj) 

The condition H18j) must hold for all Pf, and hence all j. Putting this all together yields the 
Bethe equations for the fij: 

f 

^if = i-lyl{g{^^J,^^k)■ 

k=l 

This is thus a coupled set of polynomial equations for the fii. Using the explicit form of g, and 
the expression (|19|) for the translation eigenvalue t, the Bethe equations simplify to 

^f-Vi = fT W±i^. (23) 

These are a set of / coupled polynomial equations for the Hj, j = 1 . . . f . Solving these for a set 
of Hj, one then finds the corresponding eigenstate (up to an overall normalization) by using (j22|l . 
This eigenstate has energy 

/ 

E = N-2f + Y,[H + {f^j)~'] (24) 
i=i 

The correspondence to the XXZ model at A = —1/2 discussed above is readily apparent 
in the Bethe equations (|23l) . These are precisely the Bethe equations one obtains for this XXZ 
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model with L = N — f sites, and with a twist resulting in the factor of in ()2c{|) . However, the 
eigenstates in the two cases are slightly different: the factor fJ-p{k+i)/ fJ-Pk in l|22|) does not appear 
in the XXZ case. 

In the preceding Bethe ansatz computation, we have not used the supersymmetry at all. In 
fact, one can easily repeat the computation for a model with an arbitrary coefficient of the four- 
fermion term in 1)141) . Such a model is integrable (it is related to the XXZ chain at arbitrary 
A), but it is no longer supersymmetric. It would be interesting to understand if there is a 
quantum-group symmetry generalizing the supersymmetry to arbitrary four-fermion coupling. 
In particular, it is plausible that the loop group symmetry of .9^ specialized to A = —1/2 is 
related to the supersymmetry described here. 

3.3.2 The monomer Bethe ansatz 

The Bethe ansatz has been argued to be complete [Sj- However, when / > there seem to be 
too many solutions to the Bethe equations. For example, consider the state with the maximum 
number of fermions, which has / = N/2 when is even. In this state, the hard-core fermions 
(or equivalently, the hard-core dimers) are closely packed. There are just two eigenstates for this 
value of /, but the Bethe equations are suggestive of many more solutions. Since however, all 
solutions to the Bethe equations in which the are either distinct or equal to exp(ibi7r/3) give an 
eigenvector, the number of solutions can not exceed the number of eigenvectors unless dependent 
vectors are constructed. The analogous situation in the XXZ model is easy to understand: the 
two situations are related by an overall spin flip, and the corresponding solutions to the Bethe 
equations can be related. Because of the hard-core constraint in our problem, the particle/hole 
transformation is trickier; it can be found using the mapping to the XXZ model described above, 
but will relate systems with different numbers of sites. To understand the states with / > 
it is thus useful to build up a Bethe ansatz solution starting from the close-packed state. It will 
turn out as well that the resulting Bethe ansatz equations can also be defined even in the case 
that the states have fractional charge. 

For A'^ even, there are two close-packed states. In terms of fermions they can be written as 

ici) = 44 . . . JO) \C2) = 44 • • • 4io) 

States with smaller / are described by having domain walls between the two reference states. 
A state with / = N/2 — 1 has two domain walls. It is convenient to pictorially describe each 
hard-core fermion as a dimer centered on site i. The domain walls then correspond to monomers: 
sites on the dual lattice unoccupied by a dimer. The two reference states and an example of a 
state with two monomers are then 

|Ci) = ... 

IC2) = 

mlm\Q\C) = o o 

We have denoted by rrii the operator which creates a monomer at site i on the dual lattice. 
These of course can be written in terms of the original dimer creation operators: for example, 

mimi\c) = 44444,43 ■■■m- 

All of the states in the Hilbert space can be described in this monomer basis. For even, 
there must be an even number of monomers M to ensure the periodic boundary conditions. The 
monomer construction works for odd N as well: such states correspond to odd M. For both even 
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and odd A'^, the fermion number / is / = (A^ — M)/2. Monomers can be placed next to each 
other; two adjacent monomers simply correspond to removing a dimer. The only restriction on 
monomer placement is that sequential monomers must be an odd number of sites apart. We can 
thus write any state as 

4>^^^ =(pM{h,i2,---, iM)m\^rn\^ ■ ■ ■ m\jC). 

with 1 < ii < i2 < ■ ■ ■ < iM ^ N , and ij^i — ij must be an odd number. There is no necessity 
for the subscript of \C) when there are monomers present: the monomer creation operators are 
sufficient to specify the state uniquely. 

The Hamiltonian in the monomer basis is 

A'' Ad ^ 
1=1 

A monomer can hop two sites as long as it does not hop over another monomer; this preserves 
the restriction that ij — ij+i be odd. As for dimers, the potential energy goes up when monomers 
are adjacent. 

The one subtlety in working in the monomer basis is the implementation of the translation 
operator. The translation operator takes T\Ci) = IC2), but because the dimer creation operators 
anticommute, TIC2) = — (— l)-^|Ci). In general, one picks up this minus sign if / is even and if 
the first monomer is on an even site. It is thus convenient to definite some of the monomer states 
with an extra factor of — i, namely m| m| . . . |C) = -i44 • • • |0) when both / and ii are even. 
With this definition, one has for states with no monomer on the last site A^, 

T{mlml...\C)) = Xml^,ml^,...\C) 

where A = i if / is even, and A = 1 if / is odd. If there is a monomer on the last site A^, one has 
(for any /) 

T{mlml . . . mjvIC)) = m\ml^-^ml^^ . . . m|^^_^^JC). 

The latter relation means that the monomer operators commute with each other. 
The Bethe ansatz analysis is very similar to that in the original basis. We take 

^M{ii,i2, . ■ Am) = ^^Bp /ip^/^pa • • • lipf (25) 
p 

The eigenvalue of the translation operator T is then 

M 

Imposing the periodic boundary conditions yields 

-^-SpA/,Pl,P2,...P(/-l) = ^p/-Bpi,P2,...PM- (26) 

For H25|) to describe an eigenstate of the Hamiltonian, the amplitudes Bp must obey 

2 2 1 2 

Bpi _ flpk /^PA:^P(fc+l) + 1 ~ I^P(k+l) 



Bp f^Pik+1) /^Pfc^P(fc+l) + 1 - ^Pfc 



(27) 
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where the permutation P' is defined as before. Combining ()27() with the cyclic relation ()26() yields 
the Bethe equations. In terms of new variables Vj = (/xj)^, they resemble those in the original 
basis; they are 

^^-ft = (_i)M+/-rr^<z^^fc+l^ (28) 

fell ^j^k + l-J^k 

M 

n-, = {-ly^H-^ (29) 

where we used the fact that / = (A^ — M)/2. The energy of this eigenstate is 

M 

^ = / + Eh + (^^)"'] (30) 

The similarity of the Bethe equations in the two bases is of course a consequence of the mapping 
to the XXZ model; in the XXZ model there is an obvious symmetry between up and down spins. 

3.3.3 Supersymmetry and the ground state 

In this section we describe several interesting consequences of the supersymmetry for the Bethe 
ansatz. 

The simplest one is the pairing of eigenstates of H. Acting with on the state (f>^^^ defined 
in we have 

Q+^(/) = ^ ^(i^^ . . . fj- 44A. ■ ■ ■ < 10)) . 
{i,} \i=i / 

This is a state with / + 1 fermions, and so can also be written in the form ()16() where ip takes a 
special form namely 

ii-2 12-2 N 

<^'^ = ^'P(.i,h,---'if) - ^ (p{ii,i,i2...if) ^ h (-1)-''+-^ ^ ip{ii,i2, . . . ,if,i). 

2=1 i=ji+2 i=if+2 

The minus signs are a result of the fermions anticommuting. If the state (p^^^ is of Bethe ansatz 
form H17|). the state Q'^cj)^^'^ is also of Bethe ansatz form, with 

= /ii,/i2, • • - At/)- 

The relative minus signs work out correctly from (|22|) because g{l,^j) = —1. It is easy to 
check that if a state </>*^-^^ obeys the Bethe equations (^S)) with a set Hi,H2, ■ ■ ■ ,l^f, then the set 
(1, fJ-i, fJ-2, ■ ■ ■ , ^J'f) also obeys the Bethe equations. Moreover, the two states have the same energy. 

We have thus shown that acting with the supersymmetry generator corresponds to insert- 
ing a particle with zero bare momentum, where the bare momentum pj is defined via fij = e'^^ . 
Acting with twice annihilates the state, because states with fij = = 1 for j ^ k have 
vanishing wavefunction. Thus (Q"*")^ = as required. To complete this identification of doublets 
within the Bethe ansatz one can verify explicitly that any state obeying the Bethe equations 
satisfies Q~(j)^^^ = unless one of the ^ij = 1. An E = ground state (po therefore has all fij ^ 1; 
it also has the special property that Q^cpo = even though none of the fij are 1. 
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An even deeper consequence of the supersymmetry results from the fact that the ground state 
energy is zero. We wih derive a single polynomial whose roots give the roots of the Bethe equations 
fij (or Uj) describing the ground state. Thus for the ground state we reduce a set of / coupled 
polynomial equations of order / to a single one of order /. Results very similar to those in this 
section have been derived for the XXZ chain at A = —1/2 for an odd number of sites, or for special 
twisted boundary conditions The computations of jT^ were possible because the transfer- 
matrix eigenvalues of the corresponding six-vertex model were known explicitly, through some 
clever arguments [H]. Our result thus explains why these eigenvalues could be found explicitly: 
in our formulation, the model is supersymmetric, and so the ground state energy vanishes. 

To do this computation, we need to make some definitions. We define Wj, j = 1 . . . f, via 

— a 



qwj — 1 

where q = exp(— ivr/S). In terms of w, the energy of an eigenstate is 



E = N-^ ^ (31) 



3 

Baxter's Q- function f6| is a polynomial with zeroes at w = Wk, namely 

/ 

Q[w) = - Wk). 

The Bethe equations can then be written in terms of the Wj as 

Wj-q\ ^"^ _ ^-f Q{q''^Wj) 



qwj-lj Q{q^Wj) 

Baxter defined this function for integrable lattice models because the transfer-matrix eigenvalues 
can be easily expressed in terms of it. We define an analogous function T{w) by the equation 

T{w)Q{w){l + wf'^q-^t = Q{q~^w){q~^w + if'^ - r^q~^ Q{q^w){q^w + 1)^"^ (32) 

Note that T{wk) is finite even though Q{wk) = 0, because the Bethe equations require the right- 
hand side of this equation to vanish when w = wt Qs, well. To simplify this relation, we define 
n{w) = Q{w){l + w)^-f, so that 

T{w)'R{w)q-^t = n{q-^w) - t-^q^n{q^w). (33) 

The eigenvalue t of the translation operator can be expressed as 

*-n(^^) -Q(^^ 

We have defined T{w) so that we can obtain the energy eigenvalues H24|l by expanding T{w) 
around w = q. Because TZ{q^) = TZ{—1) = 0, the relations (PS]) and give 

T{q)=t-\^^^'''"^ 



n{q) 
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Expanding T{q) and TZ{q) in Taylor series, and substituting into (jSH), relates their derivatives as 



^'C)-,-^-^ (35) 



Using the definition of TZ{w), we find that 

TZ'jw) _ N-f 1 

j=i -I 

Substituting this into (|35|) . and comparing to the energy (|31|). we find that 



T{q) g2 _ 1 



(36) 



In light of the correspondence of our model with the XXZ model discussed above, it is natural 
to expect that T[w) for arbitrary w will provide the eigenvalues for some lattice model with w 
as an anisotropy parameter. The eigenvectors of the transfer matrix of the lattice model will be 
identical to the eigenvectors of our quantum Hamiltonian. These eigenvectors will be independent 
of w, so the lattice model will have commuting transfer matrices and hence be integrable as well. 
The eigenvalues, which do depend on w, are determined by the same Bethe equations as our 
fermion chain, via the T-Q relation 1)32(1 . One can also work out a dual Tm, using the monomer 
Bethe ansatz. We will not, however, pursue this connection with the lattice model any further. 

The fact that the energy is simply related to T{q) via ()36() allows us to explicitly find Q{w) 
for the zero-energy ground states. Because = 1 and T'{q) oc E, this means that we should 
search for a Q{w) satisfying with T{w) = 1 for all w. Such a solution (denoted Qoiw)) is 
the Q- function for a ground state with E = 0. The roots of the polynomial Qoiw) are the Wj for 
this ground state. Denoting the corresponding TZ{w) by TZo{w), we require therefore 

TZoiq-^w) = tq-^TZoiw) + t-^^Uoiq'^w). (37) 

If there exists a solution to this equation with non- vanishing wavefunction, then we have an E = 
state. We will find such solutions explicitly. 

The crucial fact to exploit is that by definition Qoiw) and TZo{w) are polynomials of order / 
and N respectively. We define the coefficients so that 1Zq{w) = RfYl!k=oi~''^)^'^k, where Rf 
is an unimportant overall coefficient. Using this expansion in ()37j) gives 



for all k. Thus either the expression in parentheses vanishes, or else does. This requires that 
every third term in TZq must vanish. A non- vanishing ro means that 1 = tq-^ + t-^q^ , so for 
N = 3p with p an integer, TZq ^ only if t = (— l)^g^^, while when N = 3p ^ 1, one must have 
t = i-l)^. 

Two facts about TZo{w) give us enough information to find it exactly, following the method of 
[TO] . The first is that every third term vanishes. The second is that Qo{w) = 1Zq{w)/{1 +w)^~f 
is a polynomial of order /, so that 

iz,{-i) = Ki-i) = n'^{-i) = ... 4^-^-^)(-i) = 



14 



If / is too small, then there are too many equations to solve, and no non-trivial solution exists. If / 
is too large, one can find solutions, but these must have vanishing wavefunction (this presumably 
can be shown by utilizing the monomer Bethe ansatz). One finds a non-trivial TZq^w) only when 
f = fo = int((A^ + l)/3). Thus the ground state has / = /o, as had already been shown using 
other techniques in [H] . With a little more work, we find TZq explicitly. The non- vanishing are 

r^_3fc = (-l)^r3fc = B-^^jl'^/\k) iV = 3/0 + 1 

r7V-3fc = (-l/r-sfc = Bl^^^ff^k) iV = 3/0-1 

and 

^3fc = Bljl'l^' (k), r^k+i = Bl[^ljl {k) iV = 3/0 

where B'^'^{k) is the product of binomials 

There are two ground states for = 3/o; the other TZq is given by exchanging r^k with rN-sk- 
Dividing by (1 + w)'^~f gives Qq. There are many interesting mathematical issues associated 
with these polynomials. In particular, it is from these polynomials that some few of the very 
many observed properties of the ground state of the XXZ chain have now been proven^^l- 



3.4 Excitations 

In this subsection, we will show that the low-energy "holon" excitations of this model have 
fractional charge ±1/2. The fractional charge of 1/2 is quite natural in the monomer basis, 
because two monomers have a charge of only one. Another hint of charge-1/2 particles is the fact 
that there are two ground states when the number of sites is a multiple of three. We will use the 
Bethe ansatz to derive their existence. 

The presence of fractional charge was conjectured in by using a heuristic Neel-like picture 
of the ground state. It was suggested there that excitations with charge 1/3 appear, but there 
are several reasons why these states do not appear in the spectrum in the continuum limit. The 
reason for the charge 1/3 was that there seem to be three Neel-like ground states, differing by 
translation by one lattice site. However, the exact index results of j3| and the results for Qq 
above indicate that there are either one or two ground states, depending on the number of sites. 
Moreover, in the conformal field theory describing the continuum limit, the charge-1/3 states are 
created by a chiral operator called the spinon. Although one can build up the Hilbert space of 
the conformal field theory by acting with the spinous, neither single-spinon nor two-spinon states 
appear in the spectrum, as defined by the toroidal partition function. 

To find the charge ±1/2 states, we need to study the solutions of the Bethe equations ()23|) 
in the — > 00 limit. To make contact with past studies of the XXZ model at A = —1/2, it is 
convenient to define uj = ln(u;j)/2 so that 

_ sinh^Uj ± i7r/6) 
smh{uj — i7r/6) ' 

and the Bethe equations are 

/ sinhjuj + m/6) \ ^"^ ^ f X sinh(nj - Uk + 'nr/3) 
\smh.{uj — ivr/G) / ^--^ smh{uj — Uk — Itt/S) 
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In the N ^ oo limit, there are solutions of p8|) with Uj real, and with uj + 'm/2 real, so that both 
sides of the equation have modulus 1. There are also "string" solutions, where there are two or 
more Uj with the same real part, and having imaginary parts differing by ivr/S. The classic result 
of ^21 is that in the N ^ oo limit, the only string solution which needs to be included is the 2+ 

(2) (2) (2) 

string, which consists of pairs Uj ' — ivr/G, u - + ivr/G, with u"- ' real. Other string solutions occur 
at most in finite numbers as ^ cxd, and so can be ignored. 

Another long-known result is that the ground state of this model contains only solutions 
of ()38|) with real Uj |S] (to avoid complications involving corrections, we assume in this 
subsection that is a multiple of 3). In terms of the bare momenta pj defined by /Xj = e*^^ , real 
Uj corresponds to vr > \pj\ > vr/S, so one can think of pp = ±7r/3 as the Fermi levels. One can 
check using the following analysis that including other solutions of the Bethe equations raises the 
energy. (One can also check that all the roots of the polynomial Qo correspond to real Uj.) 

Since there are /o = int [{N + l)/3] particles in the ground state, the number of /ij in the 
ground state also diverges as ^ cxo. We can therefore define a density Pq{u), so that in this 
limit, the number of solutions of H38|) in the ground state with Uj real and in between u and 
u + du \s po{u)du. To derive po, we follow the standard procedure (see e.g. 0). We take the 
derivative of the log of (|H5|) . yielding 



po{u) = {N - fo)ai/Q{u) - / dM'po(^^')ai/3(^^ " 

J — oo 

where 

Id sinh(ii + ia) 1 sin(2a) 



u') 



27ri du sinh(u — ia) vr cosh(2M) — cos(2a) 
This can be solved by Fourier transformation, using the fact that for < a < 1/2 

^ sinh[(l - 2a)7rA;/2l 

smh[7rK/2J 

This yields 

~ ^-/o ( ^ 3(iV - /o) 



2 cosh[7rA;/6] ' 27rcosh(3n) 

This result can be checked in several ways. First, note that the fermion number /o in the ground 
state is 

/o = / dupo{u) = po(0) = — r— ^. 



Thus /o = A^/3 as shown previously. Second, one can plug in this expression for po{u) into (|24)1 . 
and verify that indeed Eq = 0. 

To understand the excitations over this ground state, we define P{u) to be the density of 
solutions of (|38() with u real. For the ground state with a multiple of 3, all real solutions are 
utilized and Po{u) = po{u). For excited states, we allow there to be holes in this distribution: a 
hole is a value Uj which solves H38() but which is not one of the u^, k = 1 . . . f. We denote the 
values of uj corresponding to these holes as u^J^^ , a = 1 . . . H. Excited states can also contain the 
two kinds of solutions of (|38j) which have u not real. Those of the form Uj = n[^^ + i7r/2 with 
u^^ real are usually called 1— strings, and are indexed 6=1.. .Mi. The 2+ strings, consisting 
of pairs Uc — ivr/G, -Uc + i7r/6, are indexed c = 1 . . .M2. The Bethe equation for the Fourier 



16 



Ke.l = /^d^xP(n) = P(0) = ^-l + ^ + 



transform P{k) for the low-energy excitations is then 

^(^) = ^°(^) + 2^;^i^(^^°^ ' -Ee j (39) 

We can use to compute the charge of the excitations. Denoting by Aureal the total number 
of real solutions of the Bethe equations, we have ()39() : 

3/0 f,H^Mi_M2 
2 2 

By definition of the total number of filled real solutions is Aureal ~ H. The fermion number / 
is therefore 

the latter factor of 2 because each 2+ string corresponds to two and hence fermion number 2. 
Thus the fermion number of a low-lying excited state is 

/ = /o-f 

One must choose the numbers of holes and strings in a given solution so that AAi-eai is an integer. 
This means there must be an even number of holes, so overall the fermion number is indeed an 
integer, as it must be. However, a two-hole state is described by two variables, u^'^^ and u^'' . 
Moreover, one can check that the energy of this state is the sum of two terms, one depending 
on u'l^^ and the other u'qI^ . Thus this state is properly described as a two-holon state with 
fermion number —1 relative to the ground state. These holons are identical, so each holon is 
interpreted as having fermion number —1/2. One can find a charge -1-1/2 state by acting with 
the supersymmetry charge Q^; this yields a state with two holes and one 1— string (the latter 
with u^^^ = i7r/2). The state has no net charge over the ground state, and is interpreted as a 
holon/antiholon pair. 

To find the excitations with positive fractional charge directly, it is easiest to look at the 
monomer basis. The analysis is virtually identical to that above. The number of monomers M 
is related to fermion number / by / = [N — M)/2. For these Bethe equations we denote Aureal, 
Mh^ A/^ii and Mq the number of real solutions, the number of holes, the number of 1— 
strings, the number of 2+ strings, and the number of monomers in the ground state. We find 

M = Aureal -MH + M1+M2 

and 

. . 3Mo M Mh Ml M2 

■^-i = ^ + T + ^ + ^ T- 

Putting these two together yields 

Mo-M Mh , 
f-fo = ^ = -^1 + M2. 

Thus we find that holes in the sea of monomers have charge 1/2. One slight difference with the 
preceding analysis is that one can find solutions of the Bethe equations where Mh = 1, so that 
the state seems to have fractional fermion number. These states {e.g. N = 15, Mq = 5, M = 4, 
A4reai = 6) do have Aureal and M integer. However, these states are not in our original space of 
states, because they violate the restriction that the number of monomers M be odd for odd N, 
and even for even N. However, it is interesting that by extending the space of allowed monomer 
states, one can find sensible Bethe equations for states with fractional charge. 
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3.5 The holons in field theory 

It is not hard to check that the low-energy excitations of the theory are gapless in the continuum 
hmit. Thus one would expect that there should be a conformally-invariant field theory describing 
the model in this limit. This conformal field theory is the simplest AA=(2,2) supersymmetric 
conformal field theory. The (2, 2) means that in this Lorentz-invariant gapless field theory, the 
symmetry is enhanced to two left-moving and two right-moving supersymmetries. Its relation 
to this lattice model has been described in detail in 0. If we assume that the charge ±1/2 
holons form a doublet under supersymmetry, this implies that the left-moving holon states are 
^±i/2,±3/4|0)a^S) ^-iid fhe right movers are V^i/2.^3/4\0)ns- The left-moving state has left and 
right scaling dimensions (3/8,0), while the right-moving state has (0,3/8). This fits in nicely 
if somewhat strangely in the M={2,2) language: the left-moving holon has left part in the first 
excited Ramond state, while the right part is in the Neveu-Schwarz vacuum. Putting a left and a 
right holon together gives the dimension (3/8,3/8) state, which corresponds to a chiral primary 
field acting on the Ramond vacua. The holon creation operators acting on the Ramond ground 
states (those appearing in the lattice model with periodic boundary conditions) are l/_j_i/2,±i/4! 
and have dimension 5/24 and Lorentz spin 1/8. It would be interesting to use the monomer 
Bethe ansatz to compute the finite-size energy of a single holon state in order to verify this. 

We note also that one can obtain a J\f=2 supersymmetric massive field theory by perturbing 
the conformal field theory by the only supersymmetry-preserving relevant operator, Vo,2 + ^,-2) 
of dimensions (2/3,2/3). This model is integrable, and its exact particle spectrum is known. In 
the N=2 language, this amounts to a doublet of charge ±1/2 kinks JJj. These kinks interpolate 
between the two vacua of the theory. It is natural to identify these two fractionally-charged kinks 
with the holons at the critical point. 

4 M2[a;]: supersymmetric models of single fermions and pairs 
4.1 The Hamiltonian and the ground states 

Our discussion in section [2.21 provides a somewhat implicit definition of the model -/Vf2[x], with 
the parameters yi, y2 fixed at yi = x, 2/2 = 1- Its Hamiltonian acts on the space 7i2 of single 
fermions (s) and nearest neighbor pairs (p) on a ID lattice, all spaced by one or more unoccupied 
sites. The Hamiltonian can be written out in second quantized form, but we prefer to list the 
non-vanishing amplitudes that it entails 

• single hop, . . . 0100 . . . . . . 0010 . . ., with amplitude — 1, 

• pair hop, . . . 01100 ... . . . 00110 . . ., with amplitude —1, 

• split-join, . . . 01010 . . . ^ . . . 01100 . . ., ... 01010 . . . ^ . . . 00110 . . ., with amplitude x, 

• partner swap, . . . 011010 010110 . . ., with amplitude 1, 

• potential, given by x^Ns + Np, where Ng counts the number of sites where, within TI2, a 
single fermion with empty neighboring sites can be put in or taken out, and, similarly, Np 
counts the number of sites where a fermion can be put in or taken out, such as to convert 
a pair to a single fermion or vice versa. 

As an example of the potential, the diagonal term for the configuration 0011010 on a 7-site 
open chain will be 2x'^ + 3. The potential can be rewritten in terms of the number of singlets A'^g, 
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the number of pairs Np, and the number of "nearest" neighbor states, where Ngs, Ngp and Npp 
counts the number of times singlets or pairs are separated by only one site in between. We have 
for N sites 

x^iV, + Np = x^N - {2x^ - 2)Ns - (4^2 - 2)Np + {x^ - 2)iV,, + (x^ - l)iV,p + x^Npp . (40) 

Clearly, x = 0, x = 1, x = \/2, and x — > oo are special in the sense that some of the amplitudes 
drop out and/or become equal to each other. 

The general discussion of section U7I\ applies to this case, so we already know the Witten index 
for any M2[x]: W2 = 3 for a chain of length N = 4n, so we know that there are at least three 
ground states. Based on numerical analysis, we expect that for general x 7^ 0, there are precisely 
three E = ground states at fermion number / = N/2. 

For X = 0, the situation is quite different and there are many more E = ground states, 
although for N = 4n there must be three more bosonic ground states than fermionic ones, so 
that the Witten index remains three. A special feature at x = is that the split-join amplitude 
in the Hamiltonian vanishes. This means that both Ng and Np are conserved individually, with 
the total fermion number equal to / = A^^ + 2Np. We observe the following ground states in 
M2[x = 0]: 

• the empty state with Np = Ng = 0, t = 1 

• at t = 1: the unweighted sum of all states with Np = 1, Ng even 

• at t = 1: the unweighted sum of all states with Np = 0, Ng odd 

• for A = 4n, at t = ±i, the states |1010 . . . 10) ± i |0101 ... 01) 

• for A = 4n + 2, at t = -1, the state |1010. . . 10) - |0101 ... 01) 

Note that the supersymmetry gives a nice way of finding the ground states: it is often easier to 
find states annihilated by both and Q~ directly, instead of those with zero eigenvalue of H. 

An unweighted sum over states amounts to having bare particles of momentum zero. Thus 
the ground states at x = include fermion condensates with any fermion number: the chemical 
potential and interactions compensate for the fermi repulsion. In the model Mi, the fermi level 
is \pf\ = Tr/3. That model has the special property that multiple fermions are allowed with the 
same bare momentum equal to the fermi level. The same property holds in M2[x = 0], but even 
more remarkably here, the fermi level here corresponds to = 0, a state with zero bare energy. 
This is why one can have an arbitrary number of such particles in the ground state, as long as 
the periodic boundary conditions are satisfied. At any finite temperature, however, there are no 
p = particles in the state minimizing the free energy: adding a p = particle to a given state 
reduces the entropy without decreasing the energy, and so raises the free energy. 

4.2 Relation with the su{2\l) symmetric tJ model when x = 

We earlier discussed a mapping between the simplest super symmetric lattice model Afi and the 
XXZ chain at A = —1/2. We can find a similar mapping for M2[x = 0]. Rather surprisingly, 
this turns out to lead to a well-known integrable model, the so-called tJ model at its integrable 
ferromagnetic point J = — 2, t = 1. 

First, it is useful to give the operators and H for x = in terms of the creation and 
annihilation operators for single fermions (s^ and s) and for the pairs {p^ and p). Expressed in 
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terms of the basic fermion annihilation operators we have Si = di = Vi-iCiVi+i at lattice site i. 
The pairs are located at the sites of the dual lattice, of which the indices take half-integer values: 
Pi+i/2 = 'Pi-iCiCi+iT'i+2- Singlets must be 2 sites from each other, pairs must be 3 sites from each 
other, and a pair and a singlet must be 5/2 sites from each other. (One can thus think of this as 
a system of hard-core dimers and trimers.) The supercharges in this basis are quite simple: they 
just change a singlet to a pair, and vice versa, namely 



N 



N 



/2 '-^1+1/2 



^1/2 



i=l 



The Hamiltonian is 



N 



H2[x = 0] = 2N, + 2Np + "^[{-sl^^Si - pl+i/2Pi-i/2 + pI+i/2P*+5/24+3«i + h.c.) 

i=l 

- '^4+2^1+24^1 ~ pl+5/2Pi+5/24si ~ ?'l-5/2?'i-5/2^1si] • (^l) 

In these formulas, the action of and H is restricted to the k = 2 Hilbert space 7i2- The 
various amplitude listed in above are in agreement with (|41j) . This model is reminiscent of the 
tJ model in that it involves two species of particles, with four-particle interactions. Note also 
that in this Hamiltonian the chemical potentials are negative and the interactions are attractive, 
both in contrast to the Hamiltonian for Mi. It thus resembles a ferromagnetic tJ model; we give 
the precise mapping next. 

The tJ model describes fermions with spin hopping on a lattice. No double occupancy is 
allowed, and the interaction between nearest-neighbor fermions is of Heisenberg type. Thus we 
have three possible states at each site, (T, i, hole). Comparing with section 1)3. 2|) . we again map 
an edge between empty sites to an up-spin, and a single occupied site to a down-spin. In addition, 
we map a pair to an empty state (hole) in the corresponding spin model. In the resulting state 
of the spin model we include a minus sign for every pair (hole at j, down-spin at k) with j < k. 
In terms of A^, Np and Ng, the parameters of the spin model are 



= N - 2Ns - 3Np, N^ = Ns, = Np 



(42) 



and the spin model lives on a chain with L = N — Ng — 2Np sites. Converting the Hamiltonian 
(|41j) to act on the spin model produces precisely the Hamiltonian of the tJ model, with J = —2 
and t = 1: 

H2[x = 0] ^ Htj[t = 1, J = -2] + 2Nh , (43) 
where the tJ Hamiltonian is given by (we follow the conventions of |2()j ) 



L 



(44) 



The operators Ql^, Qj.a and Sj and nj are bilinears of the tJ fermion annihilation/creation 
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operators Caj,clj and the bosonic hole annihilation/creation operators bj,bj at site j: 





= 4, Ad 


- ^j, 


-cr) ; Qj,(y 








Si 




Sj = ^Kt-^u 


nj,a 




Uj 




i ■ 



These operators all act within the physical space defined by rij^-^ + nj^i + bjbj = 1 for all sites j. 

As for the mapping Mi <-> XXZ, the resulting tJ model is subject to twisted boundary 
conditions, with the twist depending on the eigenvalue t of the translation operator. In particular 

-Ui,T = 4,T -Uu = (-l)^+^^-^^^^-^4a ^Ui = (-l)^^t-^4 • (45) 
The tJ model at J = ziz2t, with periodic boundary conditions, is invariant under the graded 
Lie algebra symmetry su{2\l). On a given site, the s'u(2|l) symmetry rotates the three possible 
states on the site amongst each other. The algebra is generated by Qo-, Qt, 5** and rio-. Since an 
electron is fermionic and the empty site bosonic, some of these generators are fermionic, and some 
are bosonic. For this reason graded Lie algebra symmetries are usually called 'supersymmetries' 
in the condensed-matter literature. This terminology is unfortunate, as the sii(2|l) symmetry of 
the 'super symmetric' tJ model is unrelated to the (spacetime) supersymmetry described in this 
paper. Because of the twist involved in the mapping M2[x = 0] ^ tJ, there is no Hamiltonian in 
finite volume which has both the supersymmetry and the graded Lie algebra symmetry s'u(2|l): 
in the fermion model M2[x = 0] the sii(2|l) structure is hidden, and in the tJ model there is a 
hidden supersymmetry. 

The mapping does elucidate some aspects of Al2[x = 0]. In the absence of twists in the 
boundary conditions, the tJ model with ferromagnetic sign of the spin-spin interaction has an 
E = ground state multiplet, which is a multiplet of dimension 2L -|- 1 of the su(2|l) symmetry, 
at t = 1. It consists of two su{2) multiplets, one of dimension L + 1 with highest weight state 
with N-^ = L, = Nfi = 0, the other of dimension L and highest weight state with N-^ = L — 1, 
Ni = 0, Nh = 1. In the presence of the twisted boundary conditions, eq. (|45|) . the states with 
A^l = or + odd survive as E = ground states. Mapping back to M2[x = 0] produces 
the £^ = ground states with t = 1 and Np = Ns = 0, {Np = 0, Ns odd) and (iVp = 1, even) 
that we discussed in section UlTl 

4.3 Relation with integrable spin-1 XXZ model when x = \pl 

The Hamiltonian for M2[x] for x 7^ does not conserve Ng and Np individually, but only the 
fermion number iVg -|- 2iVp. As a result, it is quite complicated. Nevertheless, at x = \/2, it is 
integrable. We show this is so in this section by mapping it onto an integrable spin-1 XXZ chain. 

This mapping proceeds similarly to the above two mappings. Because the singlets and pairs 
in M2[x] take up multiple sites, the correspondence is between M2[x] on the periodic chain with 
iV sites and the spin-1 XXZ chain with L = N — f sites and twisted boundary conditions. 
The states in the spin-1 XXZ chain are denoted +, and — , with the Hamiltonian conserving 
iV+-iV_ = L - No-2N^. 

We identify the empty state in M2[x] with the state with all spins + in the spin-1 chain. The 
singlets are identified with spin and the pairs are identified with spin — in the same fashion as 
in the tJ model above. The four off-diagonal terms in the Hamiltonian for M2[x] listed at the 
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beginning of this section then correspond respectively to the processes 



...+ 0... ^ ... + ... 

... + -... ^ ...- + ... 

...00... ^ ... + -... and ...00...^...-+... 

... - 0... ^ ... - ... 

The coefficient of these terms in the M2[x] Hamiltonian are respectively — 1, — 1, x and 1. Thus 
if we demand that the spin-chain Hamiltonian satisfy a Z2 symmetry under + <->—, this fixes 
the coupling in the corresponding model M2[a;] to be x = \pl. At this value of x the diagonal 
part of the Hamiltonian, denoted by IN^ + Np in M2[x], reduces to the simple expression 

2L - iV+_ - A^_+ 

in the corresponding spin chain, also manifestly + <-> — symmetric. This means that M2[\/2] 
indeed corresponds to a spin-1 XXZ chain. 

To show that M2[\/2] is integrable, we now show that the corresponding spin-1 XXZ chain 
is integrable. Since the 5'f7(2) symmetry of the spin-1 chain is broken to f7(l) x Z2, there 
are many possible nearest-neighbor couplings. The model is integrable only in a one-parameter 
family of these couplings These couplings can be obtained from the spin- 1/2 XXZ chain by 
implementing a procedure called fusion on the related six- vertex model with A = — cos(?7) jl5j . 
One finds that the integrable spin-1 Hamiltonian on L sites can be written as 

L 

^xxz-i = + - - - 

+ (l-2cos(r/)){5|5|+i, 5+57^1 + 57 5+,J 
+ cos(2r?)5]5]+i - cos(2ry)(5|)2(5|+i)2 + 2(cos(2,7) - 1)(5|)2 




5" = 00 5+ = 001 S 





where 



The matrices 5^ are defined without the \/2 appearing in the usual SU{2) spin-1 generators. 
When r] = 0, SU{2) invariance reappears; when rj = 7r/4, cos(2r/) vanishes and the Hamiltonian 
simplifies. 

It is now straightforward to verify that when rj = 7r/4, the above mapping takes the Hamilto- 
nian for M2[\/2] to Hxxz-i + 4L. Note again that the shift 4L is a constant in the XXZ model, 
but not in the original fermion chain. To find an equivalent model where the shift is independent 
of L, one can include a magnetic field term — Ylj=i '^^j spin-chain Hamiltonian. Then 

the Hamiltonian for M2[\/2] maps on to Hxxz-i + 2N with Tl = —2. 

The Hamiltonians of M2 [\/2] and the spin- 1 XXZ chain at A = -I/V2 are therefore identical, 
up to a shift and a twist in the latter. As with the tJ model, when one takes into account the 
periodic boundary conditions, the spin-1 model must have twisted boundary conditions. The 

diagonal terms — A'^-i — |_ clearly show that the spin model is antiferromagnetic. Despite the 

fact that this Hamiltonian is relatively complicated, implementing the Bethe ansatz is relatively 
straightforward. The fusion procedure ensures that the Bethe equations will be the same as for 
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the spin-1/2 XXZ model with A = —1/^/2] only the expression for the energy will change. The 
ground states for periodic boundary conditions will have -Eq = 0, so this means it is likely that 
So for M2[^/2] can be calculated explicitly as for Mi. 

Just as the spin-1/2 XXZ model at A = —1/2 is described by an AA = (2, 2) supersymmetric 
conformal field theory in the continuum limit, so is the integrable spin-1 XXZ model at A = 
— l/^/2. That this is the value of A which gives a supersymmetric field theory can be checked 
by studying the system at finite temperature. It is straightforward to derive the thermodynamic 
Bethe ansatz (TBA) equations giving the free energy in the continuum limit of this theory (see 
e.g. JS!)- They are identical to those for the N = (2, 2) superconformal field theory with c = 3/2, 
as can be shown by taking the massless limit of the TBA equations for the supersymmetric field 
theory with Landau-Ginzburg potential X'^ — X"^ found in |17j . 

4.4 Nested BA solution at x = 

We can solve the models M2[x = 0] and M2[x = \/2] by using the Bethe ansatz. Although this 
solution can be inferred from the mappings to the tJ model and the spin-1 XXZ model with 
twisted boundary conditions, it is instructive to solve the model directly. We do so in this section 
for M2[x = 0]. 

Since there are two different kinds of particles, the Bethe ansatz is more complicated. Particles 
of different species can hop over each other, so we need to also track of the ordering of a given 
state. We denote the locations of the singlets by ii, i2, . . . iNs-, ^-^d the pairs by zat^+i, • • • ^m, 

where M = Ng+Np is the total number of particles. A particular ordering iqi < iQ2 < ■ ■ ■ < iqm, 
is labeled by Q, a permutation of 1 . . . M. The ansatz for the wavefunction is then 

^(^^'^^) (n . . . ^M) = E Mq /^p'lVpf . • • f^X 

where P is another permutation of 1 . . . M. If all the particles are of one type, then we don't need 
the sum over Q: the symmetry or antisymmetry of the wavefunction relates the wavefunctions 
for different orderings of the particles. The eigenvalue of the translation operator is 

M 

as before. If this is an eigenstate of H, the eigenvalue is 

M 

E = 2M-^[f,j + {f,,)-']. 
i=i 

Let us first look at the system with only pairs present. We see from the Hamiltonian (|41() 
that this is a free system, except for the hard-core repulsion. The repulsion means that for any 
pair of particles, an eigenstate must satisfy 

<^(0'^'')(. . . , i + 1, i + 3, . . . ) + (^(o,7V^)(. . . , i, i + 2, . . . ) = 

for any i. This means that 

EA ( i+1 j+3 I i i+2 \ n 

p 
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for all i and j. This is satisfied if 



2 3 2 

A'p _ /^p(j+i) + f^Pjf^p(j+i) _ ^p(j+i) 
Ap + /^p(j+i)Mpj /^Pj 

This yields the Bethe equations 



Thus as before, the hard-core repulsion effectively reduces the size of the system and inserts a 
twist. 

The system with only singlets is very similar to the system Mi. Their Hilbert spaces are iden- 
tical, while the Hamiltonian for M2[x = 0] has twice the potential, and is multiplied by an overall 
minus sign. The singlet-only subsector thus corresponds to a hard-core fermion representation 
of the XXZ model at A = 1, the ferromagnetic Heisenberg model. This does not mean that the 
ferromagnetic Heisenberg model is super symmetric, because the supersymmetry converts singlets 
to pairs, and takes the theory out of the singlet-only sector. Repeating the earlier Bethe ansatz 
analysis with these slight modifications, one has 

^('''^'\...,i + l,i + 2,...) + ^('''^'\...,i,i + l,...)-2^('''''\...,i,i + 2,...) = 

These vanish if 

Dp' ^ y xy + l-2y 
Dp X xy + 1 — 2x 

where henceforth we denote fipj = x and /ip(j+i) = y for simplicity. We have denoted the 
coefficient Dp = ^p,(i2...M) here. The Bethe equations are 

The correspondence with the Heisenberg model is not surprising in light of the map to the tJ 
model; in the latter with an electron on every site, the hopping term is zero, and the model 
reduces to the Heisenberg model, with the spin up and down electrons playing the role of the 
Heisenberg spins. In the map of the M2[x = 0] to the tJ model, a state with only singlets in the 
former indeed maps to a state with an electron on every site in the latter. 

These two special cases in fact yield all the two-particle excited states, since we can use 
on ^^"^'^^ and Q~ on the to obtain any excited state with Ng = Np = 1. It is still useful 

to apply the Bethe ansatz here, to derive the i?-matrix necessary for the nesting. Let ii be the 
location of the singlet, and 12 be the location of the pair. We then write the Bethe ansatz as 



^2 



Bux'^f'' + B2iy''x''' h<i2 
Cux'^f^ +C2iy'^x'^ h>i2 



In the previous notation, Bu = ^(12), (12), -B21 = >i(2i),(i2), Cu = ^(12), (21), C21 = ^(2i),(2i)- This 
is an eigenstate of the Hamiltonian 1)41(1 if 

= (^(I'l) + + 5/2) + (^(1-1) (i, i + 3/2) - ^(I'l) (i, i + 5/2) + {i + 3,i + 1/2) 
= 5/2) + - 3/2) - - 5/2) + ip^^'^\i -3,i- 1/2) 
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These equation can be verified with the exphcit ip^^'^^ found by applying Q to (f)^^''^^ and to 
Substituting the exphcit form of ip^^'^'^ into the first of these equations yields 



= BW/\1 + xy-y) + B2ix'"\l + xy - x) + Ci2X^'^y'^ + Csiy 



while second one yields 



= Ci2y^'Hl + xy-y) + C2ix3/2(i + 



xy-x)+ Bi2X ^/^y^ + B2iy ^^^x^ 



Since there are two linear equations with four unknowns, we can solve for i?2i 5(^21 in terms of 
Bi2,Ci2, yielding 



B21 
C21 



1 



,3/2 



xy + 1 — 2x x^l'^ 



— \x 



y -^/^(x 



1) x^l'^y^l'^{x-y) 
-y) -(x-l)(y-l) 



B12 
C12 



Putting all this together gives us the i?-matrix of the theory. Denoting the vector 
{Ap, Bp, Cp, Dp), the i?-matrix is defined as .^21 = R{x, y)S,i2- Here the matrix R is 



R 



,3/2 



xy + 1 — 2x x^/2 



f^{xy + l-2x) 



V 





(x-l)(2/-l) -(xy)V2(^_y) 
-{xy)-^^'^{x -y) (x-l)(y-l) 










This will look a little more familiar if we define 

_ 6*1 + 171/2 
^ ^ ei-m/2 







y = 







.1/2 



^{xy + l-2y)J 



so 



R 



1 



,3/2 



9 + in x3/2 



ivr 



V 












ivr 

(xy)-i/2^ 




02 + i7r/2 
02 - i7r/2 



(xy)V20 
ivr 




xV2 
„l/2 






ivr 



where 9 = 9i — 02- This is the i?-matrix which gives the Boltzmann weights for a particular 
six-vertex model in an electric field; the extra factors of x and y in the last equation cancel out of 
the resulting partition function. The i?-matrix satisfies the inversion relation R{x,y)R{y,x) = 1 
as it should. 

This i?-matrix is characteristic of those arising in N=2 supersymmetric field theories, in that 
it satisfies what is known as the free-fermion condition. Denoting the six non-zero elements 
of the i?-matrix in order as a,c,b,b,c,a, it is easy to check that they satisfy ad + bb — cc = 
0. In supersymmetric field theories, the supersymmetry constrains the elements of the two- 
particle S matrix to obey this relation |17j. We have found other integrable N=2 supersymmetric 
lattice models, and their i?-matrices also satisfy the free-fermion condition. Presumably, the 
supersymmetry constrains the i?-matrix of our lattice models to be of this form, but we do not 
have a general proof of this observation. (The name free-fermion condition is a bit confusing: it 
means that a classical lattice model with Boltzmann weights given by this i?-matrix can be solved 
by using Pfaffians, without recourse to the Bethe ansatz. Such a lattice model is equivalent to a 
free-fermion theory, but our quantum theories are not: the i?-matrix appears in the nested Bethe 
ansatz, not as a Boltzmann weight.) 
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This R matrix satisfies the Yang-Baxter equation, so we can apply the nested Bethe ansatz 
to find the energy eigenstates for general {Ns,Np). Since this procedure is standard, we omit 
the details here. An excellent reference is the original paper by Yang ^Hj. Associated with each 
particle is a variable ^i, i = 1 . . . Ng + Np, so that ln(/ij)/i is its bare momentum. We then need 
the "magnon" variables i^^ to account for the extra degrees of freedom of allowing a given particle 
to be either a singlet or a pair. If we associate the i^k with singlets in a background of pairs, we 
then have k = 1 . . . Ng. The resulting nested Bethe equations are 



TT Ki^i - ^k) T-T ^ 

Ns+Nj, 

1 = TT ~ ^ (461 

f}[ fJ-iitJ-it^k + 1 - 2j^fc) 

Each solution of these equations represents an eigenstate of the Hamiltonian with the energy 

E= [2-/i,-/xri]. (47) 

1=1 

If we instead associate the magnons with the pairs so that k = \ . . . Np, the Bethe equations 

1= TT ^^(^-^^-^fe) (49) 

with E again given by (|47j) . In either set of Bethe equations, the supersymmetry doublets are 
states with identical fij and i^k, save for one extra magnon with i^^ = 0. 

These two sets of Bethe equations have both arisen in the tJ model: our second set corresponds 
to the Sutherland or 'BFF' set |191 OU] . while our first solution is the 'FBF' set first given in j2Uj . 
The two sets were shown to be equivalent in 1201 ■ Comparing the tJ equations to those for 
M2[x = 0] confirms the twist that we already described, and the relation L = N — Ng — 2Np. 
We stress that, while in the end there is agreement between the Bethe equations, the underlying 
algebraic structures are very different: our analysis here is anchored on 4 x 4 supersymmetric 
i2-matrices, while the solutions of the tJ model involve 9x9 i2-matrices associated to the graded 
Lie algebra 'u(2|l). 



5 General Mk models and the CFT connection 

We have already mentioned in sections 13.51 and 14.31 that in the continuum limit, the models Mi 
and M2[\/2] are each described by AA = (2,2) supersymmetric conformal field theories. The 
operator content of the former case is discussed in detail in ref . |S1 ; the analysis of the latter case 
is very similar. The simplest super conformal field theories (SCFTs) are called minimal models. 
They are indexed by an integer k, and have a central charge 3k /(k + 2). The field theories 
describing the continuum limit of the models Mi and M2[\/2] are the first two superconformal 
minimal models. 
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This correspondence makes it highly plausible that the models Mfe[{xj}] are also related to 
minimal models. We thus conjecture that for a specific choice of the parameters xi, . . . , the 
model Mk has a continuum limit described by A:-th M = (2, 2) superconformal minimal model. 
One piece of evidence for this conjecture is that the value of the Witten index for M^, on a closed 
chain with = n(A; + 2) is \ Wk\ = k + l, in agreement with the Witten index in the A:-th minimal 
model HI]. 

Another, very different, piece of evidence stems from the fact that there exist a finitization 
procedure which relates the Hilbert spaces of the SOFT and the lattice model. By starting 
from the chiral Hilbert space of the k-th J\f = (2, 2) SOFT, one produces finite spaces, graded 
by fermion number /, with dimensionalities identical to those of 'H): ■ The finitization scheme 
employs a quasi-particle basis of the N = (2, 2) CFT, which is related to the so-called spinon 
basis of the A:-th minimal model of su{2) invariant CFT j22j . 

We also note that the structure of the Hilbert spaces W^, involving clustering of up to k 
fundamental fermions, seems natural from the point of view of the Zk parafermions underlying 
the /c-th minimal M = (2,2) SCFT. For example, in j^S], these CFTs have been related to 
quantum Hall states with a order-A; clustering property similar to ours. 

The relation between the models Mi and is in a way similar to that between the spin 
5 = 1/2 and spin S = k/2 nearest neighbor Heisenberg models on a ID chain. In the latter 
case, there are k — 1 adjustable parameters in the Hamiltonian, which can be tuned to produce a 
critical model having the k-th minimal su{2) invariant CFT for its continuum limit. The latter 
CFTs are closely related to the k-th minimal M = (2, 2) SCFTs. 

Finally, the presence of supersymmetry enables computations like that of sect. 13.331 even at 
finite temperature. In the simplest case, at a special value of (imaginary) chemical potential, the 
Q-function for Mi in the continuum limit turns into a Bessel function |^. This simplification 
has useful consequences, because the Q-functions for a variety of integrable field theories have 
been shown to be useful in deriving the spectral determinant for the Schrodinger equation with 
certain potentials. This simplification has been shown to take place in many cases, including the 
integrable spin-A;/2 spin chains at their super symmetric points A = — cos(7r/(/c + 2)) |^. (The 
fact that these points are described by supersymmetric field theories can be shown by extending 
the arguments at the end of sect. l4.3|l . One might wonder if the simplification in general is related 
to the presence of supersymmetry indeed, the analysis of j2l] requires the supersymmetry. 
We have shown in this paper that in the spin-1/2 and spin-1 cases, the supersymmetry is present 
even on the lattice. Therefore, the fact that the analysis of |25j succeeds for any k is yet another 
indication that the models contain a point at which they are solvable and described by the 
k-th superconformal minimal model. 
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